Ambivalence of the anisotropy of the vortex lattice in an anisotropic type-II 

superconductor. 
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We present a geometry-based discussion of possible vortex configurations in the mixed state of 
anisotropic type-II superconductors, ft is shown that, if energy considerations assign six nearest 
neighbors to each vortex, two distinct modifications of the vortex lattice are possible. It is expected 
that certain conditions lead to a first order phase transition from one modification of the vortex 
lattice to the other upon varying the external magnetic field. 
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In anisotropic type-II superconductors the upper criti- 
cal field H C 2, the magnetic- field penetration depth A, and 
the coherence length £ depend on the direction with re- 
spect to crystallographic axes of the material and these 
dependencies may schematically be represented by ellip- 
soids. We consider the case where the mixed state is cre- 
ated by an external magnetic field directed along one of 
the principal axes of these ellipsoids, say the y-axis. To 
keep the discussion simple, the above mentioned quan- 
tities are assumed to be isotropic in the xy-plane and 
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C 2 . We also assume that in the fully isotropic 
case, the energetically most favorable arrangement as- 
signs six nearest neighbors to each vortex (triangular lat- 
tice). 

We use the commonly accepted notation where the 
anisotropy is expressed as 
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In this case, A^ 2 * 1 is the screening length corresponding 
to screening currents oriented along the z-axisi The 
anisotropics of A and £ in the plane perpendicular to 
the magnetic field, i.e., the xz-plane, may both be repre- 
sented by ellipses with the same eccentricities and elon- 
gated along the cc-axis, as is shown in Fig. 1. 

In magnetic fields well below H C 2, it is the magnetic 
interaction between the vortices that provides the main 
contribution to the free energy of the mixed state. This 
means that the distances between one chosen vortex line 
and its six nearest neighbors should all be equal in units 
of A, i.e., the nearest neighbors are situated along ellipses 
that reflect the anisotropy of A. As was shown in Ref. 
3, if H <C H C 2 is oriented along the xy-plane, all such 
vortex arrangements have the same free energy, indepen- 
dent of the orientation of the triangular unit cell in the 
xz-plane. It is well known, however, that part of the 
free energy arises from the interaction of Abrikosov vor- 
tices with the crystal latticed This is why we argue that 
only the two most symmetrical configurations, shown in 
Fig. 1, are of practical interest. These lattices differ by 
their orientations with respect to the x- and the z-axis. 
Although lattice 1 and lattice 2 appear as significantly 
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FIG. 1: Two possible configurations of vortices for H oriented 
along y-axis and assuming 7 = 3.5. The small ellipses rep- 
resent single vortices, their shape indicates the anisotropy of 



different, both of them correspond to the same average 
vortex density and the distances between the central vor- 
tex line (a) and its nearest neighbors (b) are also the same 
in units of A. 

The vortex lattice arrangements that follow from the 
two configurations shown in Fig. 1 are displayed in 
Fig. 2. While the vortex lattice 1 is indeed strongly 
anisotropic and its geometrical appearance is as it might 
be expected intuitively, lattice 2, on the other hand, looks 
much less anisotropic. In real space the nearest neighbors 
of the vortex line a are the vortices b and c. These vor- 
tices form a pattern with a close to triangular symmetry, 
which is emphasized by the dashed-line ellipse in the bot- 
tom panel of Fig. 2. At the same time, if we consider the 
magnetic interaction between the vortices, the distances 
must be measured in units of A and in this renormalized 
space, the nearest neighbors of vortex a are the vortices b 
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FIG. 2: Two possible vortex lattices. Although both these lat- 
tices are drawn for an anisotropy 7 = 3.5 (solid-line ellipses), 
they are equally consistent with anisotropics that differ from 
7 by a factor of 3, as is emphasized by the dashed- line ellipses 
(see text for details). 



and b', as is illustrated by the solid- line ellipse in the bot- 
tom panel of Fig. 2. Thus, although both vortex lattices 
shown in Fig. 2 correspond to the same anisotropy of A, 
the anisotropics of the corresponding vortex lattices in 
real space are quite different. The real space anisotropy 
of vortex lattice 2 may be characterized by the eccentric- 
ity e^ 2 ) of the ellipse drawn through the nearest neighbors 
of one chosen vortex line (dashed-line ellipse in Fig. 2). 
A simple calculation reveals that if for the vortex lattice 
1, e^ 1 ) = 7, the eccentricity for the lattice 2 is = 7/3. 
This means that if, e.g., the original anisotropy is 7 = 3, 
the vortex lattice 2 is perfectly isotropic in real space. 

Although the two lattices shown in Fig. 2 are quite dif- 
ferent in real space, in both cases the distances between 
the nearest neighbors in the renormalized space are the 
same in units of A. Furthermore, if Eq. (1) is exact, 
i.e., in frameworks of the applicability of the Ginzburg- 
Landau theory, the distances between the vortices are 
also equal in units of £. This makes the free energies 
of the two configurations identical not only in low mag- 
netic field limit but in the entire magnetic field range 
H c i < H < H C 2- In this case, minor contributions to 
the free energy depending, e.g., on the symmetry of the 



crystal lattice, may change the energy balance in favor 
of one or the other configuration. Our simple arguments 
cannot serve to answer the question, which of the two 
considered vortex lattices is stable. Depending on the 
particular chosen superconductor, any of the two may 
energetically be favored. 

The situation is different if the Ginzburg-Landau the- 
ory is not quantitatively applicable, for instance, at tem- 
peratures well below T c . In such cases the anisotropics 
7? = £{*v)/£(z) an d 7A = A (z )/A (xy) can be different. 
This circumstance changes the energy balance between 
the two considered vortex configurations in higher mag- 
netic fields where overlapping of vortex cores has to be 
taken into account. It is easy to check that if 7^ > 7a, 
lattice 1 has a lower free energy and vise versa. In the 
case of 75 ^ 7a, it is possible that the energy balance 
between the two vortex lattices depends on the applied 
magnetic field. If this is indeed the case, we expect a first 
order phase transition from one vortex configuration to 
the other with increasing magnetic field. 

The main result of our consideration is that the sym- 
metry of the vortex lattice does not provide unambigu- 
ous information about the anisotropy 7 of the magnetic 
field penetration depth. It turns out that each vortex 
configuration is consistent with two unequal values of 7, 
differing by a factor of 3. For instance, vortex lattice 1 
equally well corresponds to 7 = 3.5 and 7 = 10.5, while 
lattice 2 reflects the anisotropies of 3.5 and 1.17, as it is 
shown by the solid and dashed-line ellipses in Fig. 2. 

Recognizing the fact that the anisotropy of the vortex 
lattice in real space may be quite different from that of A 
is rather important for the correct interpretation of ex- 
perimental observations. With this in mind we recall the 
recent experimental observation of the vortex lattice in 
MgB2, invoked by magnetic fields oriented along the ab- 
planes of the hexagonal crystal latticed A hexagonal vor- 
tex lattice corresponding to an eccentricity e = 1.19 was 
observed in these experiments. This value of e is much 
smaller than the corresponding anisotropy 7 of H C 2- Al- 
though some theoretical justifications for the anisotropy 
of the magnetic field penetration depth 7a = \^ /\^ xy ^ 

to be smaller than 7ff = h£ v) /H$ = were 
offered in Ref. 0, we argue that the observation of Ref. 
most likely reflects a A-anisotropy 7a = 3e « 3.6 (see 
the bottom panel of Fig. 2). This latter value is in fair 
agreement with the observed anisotropy of H C 2 £ 

In conclusion, we showed that for each value of the 
anisotropy of the magnetic field penetration depth 7a 
there are two possible arrangements of the hexagonal vor- 
tex lattice with corresponding eccentricities e' 1 ' = 7a and 
e (2) = 7 A /3. This is why an unambiguous determination 
of 7a from experimental observations of the vortex lattice 
is only possible if an approximate value of 7a is a priori 
known. 
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